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^^ ' gas with recent LGT results at finite chemical potential. We argue that for T < 

■^ , Tc the equation of state derived from Monte-Carlo simulations of two quark- 

^^ ' flavor QCD at finite chemical potential is consistent with that of a hadron 

fH , resonance gas when applying the same set of approximations as used in LGT 

Q_i' calculations. We indicate the relation of chemical freezeout conditions obtained 

' , from a detailed analysis of particle production in heavy ion collisions with the 

^j' critical conditions required for deconfinement. We argue that the position of 

r^ , a hadron-quark gluon boundary line in temperature chemical potential plane 

can be determined in terms of the resonance gas model by the condition of 
fixed energy density. 
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1. Introduction 

A detailed analysis of particle production in heavy ion collisions has shown that in a 
broad energy range from SIS through AGS, SPS^ up to RHIC particle yields resem- 
ble that of chemical equilibrium population with respect to the partition function 
of hadron resonance gas (HRG) [Q. 

In Fig. 1^ we show an example of the comparison of HRG-model with ex- 
perimental data on different particle yields ratios measured in Au-Au collisions at 
RHIC energies. The model results [ |3| E] are shown as horizontal lines whereas 

^The most recent thermal analysis of SPS data was presented in Ref. f[2i. 
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experimental data are indicated as points. The apparent quite good agreement of 
the model predictions with the experiment indicates that production of particles in 
heavy ion collisions is of thermal origin [ ^ . A collision fireball at chemical freeze- 
out is characterized by only two thermal parameters, the temperature T and baryon 
chemical potential^ /is. The volume of the fireball can be possibly connected with 
the HBT source size [ISl 171. 

Figure Q) shows the compilation of the chemical freeze-out parameters that 
are required to reproduce the measured particle yields in central Au-Au or Pb- 
Pb collisions at SIS, AGS, SPS and RHIC energy. The GSI/SIS results have the 
lowest freeze-out temperature and the highest baryon chemical potential. As the 
beam energy increases a clear shift towards higher T and lower ^lb occurs. There 
is a common feature to all these points, namely that the average energy {E) per 
average number of hadrons (TV) is approximately 1 GeV. A chemical freeze-out in 
A-A collisions is thus reached [^ when the energy per particle {E)/{N) drops below 
1 GeV at all collision energies. The physical origin of the above freezeout condition 
requires dynamical justification. Recently, this point has been investigated in central 
Pb-Pb collisions at the SPS in terms of the Ultra-relativistic Quantum Molecular 
Dynamics model (UrQMD) [[HI EI]- A detailed study has shown that there is 
a clear correlation between the chemical break-up in terms of inelastic scattering 
rates and the rapid decrease in energy per particle. If (E) /{N) approaches the value 
of 1 GeV the inelastic scattering rates drop substantially and further evolution is 
due to elastic and pseudo-elastic collisions that preserved the chemical composition 
of the collision fireball. Following these UrQMD results one could consider the 
phenomenological chemical freeze-out of {E)/{N) ~ 1 GeV as the condition of 
inelasticity in heavy ion collisions. 
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Fig. 1. Comparison of the experimental data on different particle multiplicity 



ratios obtained at RHIC at ^/s 



NN 



130 and 200 GeV with thermal model calcu- 



lations. The thermal model analysis is from Refs. [El 0]. 



^The quark chemical potential ^q = ■:7fJ,B 
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Fig. 2. Chemical freezeout in heavy ion collisions. The full line represents the 
condition of fixed energy/particle ~ 1.0 GeV from Ref. [|H|. The broken line 
describes the chemical freeze-out conditions of fixed total density of baryons plus 
antibaryons, nb + ni = 0.12/fm^ from Ref. [E]. The points are the average values 
of Tf and ^^ obtained from the analysis of particle yields in heavy ion collisions at 
the SIS, AGS, SPS and RHIC energy [p. 

Chemical freeze-out in heavy ion collisions was recently proposed to be deter- 
mined [ ^5 by the condition of fixed density of the total number of baryons plus 
antibaryons F = {ub +'^_b) — 0-12/ fm^. It is seen in Fig. |(2J) that, within statisti- 
cal uncertainties on the freeze-out parameters the above condition provides a good 
description of experimental data from the top AGS up to RHIC energy. However, 
at the low SIS energy it overestimates the freeze-out temperature for a given chem- 
ical potential. Consequently, e.g. the yield of pion/proton or pion/participant and 
strange/non-strange particle ratios obtained at SIS turn out to be larger than the 
experimental values. Taking the value of F ~ 0.05/ fm^ at SIS energy and using 
the freezeout condition of Fsis ~ fixed leads eg. at RHIC to the freezeout temper- 
ature of ~150 MeV, the value which is lower than obtained from the particle yield 
ratios as seen in Fig. ^. At fiB — the {E)/{N) condition converge to Tf ~ 170 
MeV which is well consistent with the critical temperature T = 167(13) [ 12 and 
T = 173(8) [13 expected in (2+1) and 2-flavour QCD respectively. 

At SPS and RHIC the freeze-out parameters, the temperature Tf and the en- 
ergy density e/, predicted by the presently used thermal model of hadron resonance 
gas, agree well with the values of the corresponding parameters required for de- 
confinement in LGT calculations [^ ^M ^j- The above quantitative agreement 
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of freeze-out and critical parameters suggests that at SPS and RHIC the chemical 
freeze-out appears in the near vicinity of the phase boundary [I18|. If this is indeed 
the case then the phenomenological statistical operator of hadron resonance gas 
Zhg should also provide, consistent with LGT, description of QCD thermodynam- 
ics in confined, hadronic phase of QCD ["151 11^)1 1171 119| . In the following we argue 
that the basic qualitative properties of Zhg resulting from its dependence on T 
and baryon chemical potential ^b are indeed present in the recent lattice results 
of two flavor QCD at the finite chemical potential. We also show that imposing 
a fixed energy density condition for deconfinement the hadron resonance gas par- 
tition function Zhg can also describe [ I15| the quark mass and number of flavor 
dependence of the lattice critical temperature at jjlb = and the position of the 
phase boundary in the (T, /iB)-plane at small /is [|2()j. 
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Fig. 3. Left-hand figure: temperature dependence of AP/{xqT'^) ratio for two 
different values of fJ.q/T. The 2-flavor lattice results are from [ EHl ^M- The 
lines are the hadron resonance gas model values with (broken-line) and without 
(fulHine) the Taylor expansion of cosh(/.iB/T)) [El- The right-hand side shows 
the temperature dependence of the ratios of the second, fourth and sixth order 
coefficients in the Taylor expansion of thermodynamic pressure. The points are 
lattice results from [|2S1 ^M- The lines at T < Tq are the values obtained from the 
expansion of cosh(/iB/r)). 



2. QCD thermodynamics of hadronic phase 



The phenomenological partition function used in the description of particle pro- 
ductions in heavy ion collisions was, following Hagedorn, constructed as a non- 
interacting hadronic gas which is composed of all hadrons and resonances. In this 
approach the interactions between hadrons are included to the extent that the ther- 
modynamics of an interacting system of elementary hadrons is effectively approxi- 
mated by that of mixture of ideal gases of stable particles and resonances [|21 121] • 
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In the Boltzmann approximation, suitable for the moderate values of fis < "t-jv 
and T > 50 MeV, there is a factorization of T and hb/T dependence in relevant 
observables characterizing a baryonic sector of the system. The basic quantity is 
the pressure AP = P[T, ^b) - P{T, hb = 0) 

^^F(T)(cosh(^-l) with F{T)^ j dmp{m)C^fK,{^) (1) 

from which the baryon density ns and the baryon number susceptibility A^s are 
obtained as the first and second order derivatives with respect to /is, respectively. 
In the phenomenological application of the above model the mass spectrum p{ni) 
in Eq. |^ is approximated as the sum of delta functions. The baryonic pressure 
gets the contributions from all known baryons and resonances. 

The net baryonic pressure AP{T, fig) has been recently obtained on the lattice 
in two flavor QCD as the Taylor series in Pq/T 

up to 0(/i^) order [ESI ^. 

The basic observables characterizing baryonic contribution to the QCD thermo- 
dynamics, the net quark density Uq and susceptibility Xq are obtained from Eq.lPJ 
as 

X^ ^ d^Ap{T,f,q)/T^ 

T2 d {pq/Tf 



2c2(r) + 12c4(T)(^) +30C6(T)(^) . (3) 



The coefficients ct (T) in the Taylor expansion (|2I3|I were calculated through the 
Monte-Carlo simulations of 2-flavor QCD [ 23 , 24 . 

The lattice results (|2I3|I restricted to the confined phase of QCD can be directly 
compared with the predictions of the hadron resonance gas model ^. This, in 
general, required a numerical analysis. However, there are particular features of the 
resonance gas partition function that could be directly verified by the LGT results. 

First, as seen from Eq. |^ there is a factorization of T and [pq/T) dependence 
in the thermodynamic pressure and associated baryonic observables Uq and Xq- 
Thus, in the resonance gas model any ratios of these observables, calculated at fixed 
quark chemical potential {p,q/T) = fixed, should be independent of temperature [ 
1151 I19| . In Fig. (3-left) we show as an example the ratio of AP/A{xbT'^) for two 
different values of /Xg/T as the function of T. It is clear from Fig. (3-left) that the 
factorization predicted by HRG partition function Zhg is also seen in LGT results 
for T < Tc. 

Second, the (/iq/T)-dependence in the HRG-model is described by the (cosh)- 
function. Thus, performing the Taylor expansion of the pressure 1^ in (pq/T) 
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^^i^-^'^nnCim'^'-^m'^'^m'). (4) 



T4 ^ ' \2 \T J 4! Vr/ 6! VT 

one gets the predictions on the values of the coefhcients Ci in Q from the HRG- 
model 



C2{T) = ^F{T) , C4 = ^F(r) , C6 = -F(T). (5) 

where for a discrete hadronic mass spectrum 

^(r) = E^(f)'^^(f)- (6) 

with the sum taken over all baryons and baryonic resonances of mass rrii and the 
spin-isospin degeneracy factor di. 

In Fig. (3-right) we show the ratios C4/C2 and cq/c^ obtained on the lattice 
in 2-flavour QCD and the corresponding results of the HRG-model from Eq. JSJ: 
C4/C2 = 3/4 and cq/c4 = 0.3. The temperature dependence of the Taylor coefficients 
in Eq. Q and is controlled by a common function for all Cj. Thus, in HRG- 
model, the ratios of different Ci are independent of T. The lattice results for T < T^, 
are, within statistical errors, consistent with this prediction. The values of ce/ci 
and C4/C2 ratios obtained recently on the lattice [ I23L ^^ are also seen in Fig. (3- 
right) to coincide with HRG-model. At temperature T ~ Tc, however, the lattice 
results deviate from the resonance gas values. This is to be expected as at the 
critical temperature, deconfinement releases the color degrees of freedom^ which 
are obviously not present in the statistical operator of the hadron resonance gas. 
The partition function (^ is not sensitive to the phase transition thus, could be at 
most applicable below deconfinement. 

The pressure calculated on the lattice (see Fig. (4) ) increases abruptly when 
approaching deconfinement transition from the hadronic side. If the phenomeno- 
logical statistical operator Zhg is of physical significance then this increase could 
be due to the resonance formation. To check the importance of resonances near 
deconfinement one would need to reproduce lattice results on the T-dependence 
of AP at fixed Hq/T. However, to quantify this dependence one needs to imple- 
ment the same set of approximations in Eq.(l) as those being used on the lattice. 
The lattice results were obtained in 2-flavor QCD, thus there is no contribution of 
strange baryons in Eq.(l). In addition due to the Taylor expansion of AP in the 
lattice calculations one also needs to apply a similar approximation in Eq.(l). 

^ These color degrees of freedom could be also possibly described as resonance bound states [ 
1^ . 
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Fig. 4. Temperature dependence of baryonic pressure for different values of iiqjT . 
The lattice results are from [|2S|- The lines are the statistical model results [I15|. 

The current lattice results in Fig. (4) were obtained with a large quark mass 
corresponding to to^ ~ 770McV. This large quark mass distorts the baryon mass 
spectrum in Eq. (^. The pion mass dependence of baryon masses can be obtained 
from lattice calculations at zero temperature [E| I^HI- We use the following ansatz 
for the parametrization of baryon mass dependence on the pion mass, 



TO*(?n7r) 



1 



TO 






(7) 



with A ^ 0.9 ±0.1, m* being a distorted hadron mass at fixed to^ and to is its 
corresponding physical value. 

In Fig. (4) the lattice results are compared with Eq. (^ calculated with modified 
baryon mass spectrum following Eq. Q. The T dependence of QCD thermody- 
namics obtained on the lattice at different values of fJ.q/T is seen in Fig. (4) to be 
well consistent with the predictions of a HRG~model. A similar agreement was also 
found for other relevant thermodynamic observables in (2+1), 3-flavour [^] and 
4-flavour QCD [El- The above agreements indicate that the phenomenological 
partition function of hadron resonance gas should be a good approximation of the 
QCD statistical operator in the hadronic phase. 



3. Critical conditions for deconfinement 

The recent lattice results on QCD thermodynamics show that deconfinement occurs 
at similar values of the energy density in QCD with light quarks as well as in the 
pure gauge theory, although the transition temperature shifts by about 40% and 
e/T^ differs by an order of magnitude [ E]- Thus, it has been suggested that 
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Fig. 5. The transition temperature vs. pion mass obtained in lattice calculations 
and lines of constant energy density calculated in a resonance gas model [ .iBa . The 
left hand figure shows a comparison of constant energy density lines at 1.2 (upper), 
0.8(middle) and 0.4(lower) GeV/fm'^ with lattice results for 2-flavor QCD obtained 
with improved staggered [ E] as well as improved Wilson [ I27L I28L 1^ fermion 
formulations. Tc as well as mps are expressed in terms of the corresponding vector 
meson mass. The right hand figure shows results for 2 and 3 flavor QCD compared 
to lines of constant energy density of 0.6 GeV/fm'^. Here Tc and mps are expressed 
in units of ^/a. For a detailed description see Ref. [ Q^ . 

for arbitrary quark masses the transition in QCD appears at the roughly constant 
energy density. Such an assumption can be verified by the HRG-model formulated 
for an arbitrary value of the quark mass [ fl5\. In Fig. (5) we show the lines of 
the constant energy density calculated in the resonance gas model and compare 
these to the transition temperatures obtained in lattice calculations. As can be 
seen the agreement is quite good up to masses, m^r — 3 ^/a or ttItt ~ 1.2 GeV. At 
larger values of the quark mass the observed deviations are due to the fact that 
thermodynamics is dominated by glueballs which are not present in the hadronic 
mass spectrum [ 15 . 

The critical conditions of fixed energy density could be also verified at the finite 
chemical potential.^ Fig. (6) shows recent lattice results on the position of the phase 
boundary line in 2-fiavor QCD obtained within the Taylor approximation and with 
the pion mass TOtt — 770 MeV. The line of fixed energy density e ~ 0.6 GcV / fm? 
with e calculated in the HRG-model with the modified mass spectrum is seen to 

^The critical condition of fixed energy density can be related with the percolation condition for 
deconfinement f 1301 . 
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Fig. 6. Lattice results [123 on the phase boundary curve (line with errors) together 
with phenomenological freeze-out values of T and ^b (points) obtained from the 
analysis of particle production in heavy ion collisions [T . Short-dashed and dashed- 
dotted lines are the statistical model results obtained under the condition of fixed 
e ~ 0.6 GeV/fni^ with m^ ~ 0.77 GeV and m^r — 0.14 GeV respectively. Also 
shown (full-line) is the phenomenological freeze-out curve of fixed energy/particle~ 
1 GeV from [iSj. 



coincide with lattice results. Decreasing the pion mass to its physical value and 
including a complete set of resonances expected in (2-|-l)-flavor QCD results in the 
shift of the position of the phase boundary line towards a phenomenological freeze- 
out condition of the fixed energy/particle~ IGeV. The splitting of freeze-out and 
phase boundary line appears when the ratio of meson/baryon multiplicities reaches 
the unity. 

The dashed-dotted line in Fig. (6) could be a physical boundary line separating 
the hadronic from quark-gluon plasma phase. However, in the baryon dominating 
medium one expects a very efficient meson baryon scatterings [I31|. Consequently, 
the particle dispersion relations and spectral functions can be modified in a medium. 
Thus, it is not clear if in this case the partition function of hadron resonance gas 
with free particle dispersion relations is still an adequate approximation of the QCD 
thermodynamics just before deconfinement. In view of the above, the position of 
the critical curve shown in Fig. (6) could be shifted in the parameter range where 
fiB > 0.5 GeV, that is where baryons start to be dominating degrees of freedom in 
a thermal fireball. 
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4. Conclusions 

We have shown that the statistical operator of hadron resonance gas used to de- 
scribe particle yields in heavy ion collisions provides a satisfactory description of 
the recent lattice results on QCD thermodynamics at the finite chemical potential. 
In particular, the basic property of this operator like e.g. a factorization of tem- 
perature and chemical potential dependence is obviously confirmed by the lattice 
results. The ratios of the coefficients in the Taylor expansion of a thermodynamic 
pressure obtained in LGT are as well described by the expansion of cosh-function 
predicted by the hadron resonance gas. These results should be independent from 
the particular choice of the quark mass in lattice calculations and are also to a large 
extent free from lattice artifacts. The phenomcnological partition function, with 
m^-dependent hadronic mass spectrum, was also shown to describe quantitatively 
the temperature and chemical potential dependence of the basic thermodynamical 
observables obtained in LGT for T < Tc- This indicates that the hadron resonance 
gas partition function is a good approximation of QCD thermodynamics in the 
hadronic phase. Applying the above partition function together with the condition 
of a fixed energy density we have discussed a possible position of the QCD critical 
curve in the temperature-chemical potential plane. 

Acknowledgments 

Stimulating discussions with P. Braun-Munzinger, S. Ejiri, F. Karsch, H. Satz, J. 
Stachel and A. Tawfik arc kindly acknowledged. This work was partly supported 
by the KBN under grant 2P03 (06925) and by the DFG under grant GRK 881. 

References 

1. For review see: P. Braun-Munzinger, K. Redlich and J. Stachel, in Quark 
Gluon Plasma 3 (Edts. R. Hwa and X.-N. Wang); nud-th/0304013^ 

A. Andronic and P. Braun-Munzinger, hep-ph/0402291 

2. J. Manninen, F. Becattini, A. Keranen, M. Gazdzicki and R. Stock, 
|nucl-th/0405015 

3. P. Braun-Munzinger, D. Magestro, K. Redlich, and J. Stachel, Phys. Lett. 
B518 (2001) 41 and references therein. 

4. D. Magestro, J. Phys. G28 (2002) 1745. ^_^ 

5. See also, M. Gyulassy and L. McLerran, nucl-th/0405013 



6. G. Agakichiev et al., CERES/NA45 Collaboration, 
Phys. Rev. Lett. 92 (2004) 032301. 

7. H. Satz, hep-ph/0405051 

8. J. Cleymans, and K. RedUch, Phys. Rev. Lett. 81 (1998) 5284. 

9. M. Bleicher and J. AicheUn, Phys. Let. B530 (2002) 81. 



From chemical freezeout to critical conditions in heavy ion collisions 11 

10. see also: L. Bravina, et al., Nucl. Phys. A698 (2002) 383; Phys. Rev. C66 
(2002) 014906. 

11. P. Braun-Munzinger and J. Stachel, J. Phys. G28 1971 (2002). 



12. C. Bernard et al., MILC Collaboration, |hep-lat/0405029| 

13. F. Karsch, E. Lacrmann, A. Peikert, Nucl. Phys. B605 (2001) 579; F. Karsch, 

E. Laermann, A. Pcikcrt, Phys. Lett. B478 (2000) 447; F. Karsch, Nucl. 
Phys. A698 (2002) 199c. 

14. J. Stachel, Nucl. Phys. A654 (1999) 119c. 

15. F. Karsch, K. Redlich and A. Tawfik, Eur. Phys. J. C29 (2003) 549; 

F. Karsch, K. Redlich and A. Tawfik, Phys. Lett. B571 ( 2003) 67. 

16. M. D'Eha and Maria-Paola Lombardo, GEF-TH-2004-07; "hep-lat/ 04060 12| 

17. M. D'Elia, Maria-Paola Lombardo, Phys. Rev. D67 (2003) 014505; 



Maria-Paola Lombardo, hep-lat/0210040 



18. P. Braun-Munzinger, and J. Stachel, Nucl. Phys. A606 (1996) 320; Nucl. 
Phys. A683 (1998) 3. 

19. S. Ejiri, et al., to appear. 

20. F. Karsch, K. Redlich and A. Tawfik, |nuci-th/0404009'| A. Tawfik, et al., to 
appear. 

21. R. Hagedorn, Nuovo Cimento 35 (1965) 395; R. Hagedorn, Thermodynamics 
of strong interactions, CERN Report 71-12 (1971). 

22. R. Vanugopalan and M. Prakash, Nucl. Phys. A546 718 (1992); V. Koch, 
Nucl. Phys. A715 108 (2003). 

23. C. R. AUton, S. Ejiri, S. J. Hands, O. Kaczmarek, F. Karsch, E. Laermann 
and C. Schmidt, Phys. Rev. D68 (2003) 014507; see also, S. Ejiri, poster 



session at Quark Matter 04; S. Ejiri, |hep-lat/0401012 

24. C. AUton, et al., to appear. 

25. E. Shuryak, hep-ph/0405066 and references therein. 

26. M. Gockeler, R. Horsley, D. Pleiter, P. Rakow and G. Schierholz Phys. Lett. 
B532 (2002) 63; C. M. Maynard, Nucl. Phys. Proc. Suppl. 109 89. 

27. C. Bernard, et al., Phys. Rev. D56 (1997) 5584, and references therein. 

28. A. AH Khan et al. (CP-PACS), Phys. Rev. D63 (2001) 034502. 

29. R.G. Edwards and U.M. Heller, Phys. Lett. B462 (1999) 132. 

30. V. Magas and H. Satz, Eur. Phys. J. C32 (2003) 115. 

31. W. Weinhold, B. Friman and W. Norenberg, Phys. L ett. B433, (1998) 236; 
M. Post, S. Leupold and U. Mosel, "nucl-th/0309085' M.F.M. Lutz, G. Wolf 
and B. Friman, Acta Phys. Hung. New Ser. Heavy Ion Phys. 17 (2003) 313; 
Prog. Theor. Phys. Suppl. 149 (2003) 152; G.E. Brown, Guo-Qiang Li, R. 
Rapp, M. Rho and J. Wambach, Acta Phys. Polon. B29 (1998) 2309; G. E. 
Brown and M. Rho, Phys. Rept. 269 333 (1996). 



